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ABSTRACT 

We propose that the presence of additional planets in extrasolar planetary systems 
can be detected by long-term transit timing studies. If a transiting planet is on an 
eccentric orbit then the presence of another planet causes a secular advance of the 
transiting planet's pericenter over and above the effect of general relativity. Although 
this secular effect is impractical to detect over a small number of orbits, it causes long- 
term differences in when future transits occur, much like the long-term decay observed 
in pulsars. Measuring this transit-timing delay would thus allow the detection of either 
one or more additional planets in the system or the first measurements of non-zero 
oblateness (J2) of the central stars. 
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1 INTRODUCTION 

The study of long-term orbital precession was one of the 
triumphs of celestial mechanics, when the planetary theo- 
ries of Laplace and Lagrange showed that essentially all the 
known long-term precessions of the planetary orbits could 
be explained by their mutual gravitational interaction. The 
perturbation caused by the small planetary masses 'breaks' 
the perfect central force character of the Sun's gravitational 
held, causing the planetary orbital nodes to regress and their 
perihelia to slowly advance, with typical periods of 10 4 -10 5 
years. It was this advanced understanding of celestial me- 
chanics that permitted LeVerrier and Adams to detect a 
new planet in our Solar System (Neptune) by inverting its 
observable effect on the known planets to predict Neptune's 
mass and position. 

The only known exception to these predictions at the 
start of the 20th century was the mystery that Mercury's 
perihelion longitude advanced 0.42"/yr faster than the pre- 
dicted rate of ~5.3l"/yr produced by the perturbative ef- 
fects of all the other planets. The obvious possibility was 
that the Sun's mass distribution was not spherically sym- 
metric, but explaining the mercurian advance would require 
a solar oblateness (measured by the parameter J2) of sev - 
eral percent, which was uncomfortably large (e.g. lHalll900l ). 
General relativity provided the solution, explaining essen- 
tially the entire discrepancy. In fact the Sun's oblateness is 
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left with only an empirical upper bound of J2 < 1 x 10 , 
with a current theoretical estimate for its true value of of 
J 2 ~ 1 x 10" 7 jPireaux fc Rozelodlioolh 

As of the beginning of 2006, nearly two hundred extraso- 
lar planets have been disc overed and several o f them exhibit 
transits (as catalo gued by [Butler et al.ll2006l and refer ences 
therein, especially iMavor et al.ll2004 IVogt et alj|2005j ). The 
longitude of periastron of eccentric hot Jupiters (3-day peri- 
ods) will show a secular advance of the instant of transit. If 
the orbital period can be well established (via transit tim- 
ing or radial velocity) then the long-term drift of the transit 
centers will allow one to measure the slow advance of pe- 
riastron. Observations on 10-year baselines should certainly 
show the relativistic advance, which is much more important 
for close-in hot Jupiters since they are closer to their parent 
stars than Mercury is to our Sun. In this paper we discuss 
the possibility of using the periastron advance rate to mea- 
sure either a host-star oblateness or the presence (and thus 
discovery) of additional planets in the system. 

IMiralda-Escudel (|2002h provided the hrst estimates of 
the various cont ributions to the p eriast ron advance of ex- 
trasol ar planets. lAgol et"ail (J2005) and iHolman fc Murravl 
|2005l ) build upon this work by including the effects of reso- 
nances and including other contributions to the timing noise 
of planetary transits. Although these numerical studies au- 
tomatically include the secular advance of periastron, they 
focus more on the stoci astic variation of t he int erval between 
transits. This work as IMiralda-Escudel (|2002f ) did focusses 
specifically on the secular periastron advance and builds 
upon that work by including a more accurate calculation of 
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the advance for planets whose orbits have similar semima- 
jor axes and by outlining several techniques to measure the 
periastron advance and their associated precision (including 
the effects of general relativity) . 

In 321 we calculate the secular advance of periastron of 
the orbit of a planet around a star due to other planets in the 
system under reasonable approximations. We also present 
the relativistic and quadrupole contributions to the perias- 
tron advance. ^3] places these calculations in the context of 
extrasolar planetary systems. Sj4] presents analytic and nu- 
merical estimates of how well we can determine the advance 
of periastron using various techniques (timing of the primary 
and secondary transits with or without radial velocity infor- 
mation). The special relativistic corrections are outlined in 
Sj5] and Sj6] outlines the prospects of this technique using the 
planets discovered so far as a guide. 



2 SECULAR ADVANCE 

In the Newtonian two-body problem the Laplace Runge- 
Lenz vector (also known as Hamilton's vector, or the eccen- 
tricity vector), which points from the star to the planetary 
orbit's pericenter, is stationary, so the location of periastron 
is constant. Several effects can cause the periastron to ad- 
vance. In our Solar system in increasing order of importance 
we have : 

(i) Stellar oblateness, 

(ii) General relativity, and 

(iii) Other planets. 

The amplitude of these effects depends on several parame- 
ters, and thus in extrasolar planetary systems the order of 
importance may differ. In the following sections we provide 
expressions for the rate of periastron advance for these ef- 
fects. 



2.1 Stellar contributions 

The central star can cause periastron advance by either 
being non-spherical or due to general relativistic effects 
caused by its mass. These effects cause a periastron advance 
|Misner et al.lll973h of 



and for the stellar oblateness 
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per radial period due to the star itself. Here, c is the speed 
of light, G Newton's gravitational constant, a the semimajor 
axis of the planet's orbit, e its eccentricity, and M*, R t , and 
Ji are the mass, radius, and oblateness parameter of the 
star, respectively. 
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change for the relativistic component of 
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where we have scaled these effects to typical values appro- 
priate for a hot Jupiter. For Mercury the relativistic effect 
is 0.43"/yr. This drops as a -5 ^ 2 for more distant orbits, be- 
ing 0.086"/yr for Venus, 0.038" /yr for Earth, and 0.013"/yr 
for Mars. The scaling value of J 2 = 10 -6 used in Eq. [3] is 
essentially a firm upper limit for the oblateness of our sun 
that can be obtain from theoretical and observational con- 
siderations (|Pireaux fc Rozelot|[2fJ03l ). 



2.2 Planetary contributions 

Just as Mercury's pericenter advance is affected by the other 
planets, an exoplanet's orbit will precess due to the pertur- 
bations of an unseen planet. We will generally assume that 
a second planet in the system is an external one, although 
the theory is almost identical if the perturber is interior. 

A simple way to estimate the contribution to the perias- 
tron advance from another planet in the system is to assume 
that the mass of the second planet is smeared out over a cir- 
cular ring coplanar with the first planet's orbit and calculate 
the contribution of this ring to the potential. If the orbit of 
the second planet is elliptical or not coplanar with the ob- 
served planet, the eccentricity, inclination and longitude of 
the ascending node will also change on a secular timescale. 

An additional planet in the system causes the periastron 
of an observed planet to advance by 
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during each radial orbit (|Price fc R ush 1979). The approx- 
imation holds if a <C a 2 and a 2 is the radius of the second 
planet's orbit. The mass of second planet is 7712. 

We can do a bit better if we relax the assumption of 
a <C a 2 . In this case we find that the potential due to the 
second planet is 



V 2 = 



Gm 2 



K 
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where K is the com plete elliptic in t egral of the first kind 
and a = a/a2. From IPrice fc Rush! ljl979h , the advance of 
periastron per radial orbit is 



2ip - 2tv 



(6) 



where the increase in azimuthal angle for half a radial orbit 



ij = n{3 + a[V"(a)/V'(a)]} 



-1/2 



(7) 



and V(a) is the potential due to the star and the perturbing 
planet. The limit V(a) — 1/a for no perturbing planet yields 
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ip = 7r as expected. For mi <C M* we have 
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where i? is the complete elliptic integral of the second kind. 
The function &P/ 2 (a) is a Laplace coefficient from classical 
perturbation theory; it is equal to 3a for a <C 1, is order 
unity for q=0. 1-0.5, and then increases rapidly to >100 for 
a > 0.9. In the limit of a large ratio between the two semi- 
major axes we obtain 
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a< 1, exterior perturber 
«>1, interior perturber 
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This is perhaps more transparently expressed by noticing 
that N p — 2n/ Szu is the number of inner planet revolutions 



required for a full precession of its orbit 
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for ai < d2. Fig. [T] gives the exact value of N p x (mi/M,) 
and its asymptotic form. As examples, an unseen Jupiter- 
mass planet (rni/M* ~ 10 _ ) at ten times the semimajor 
axis of an interior hot- Jupiter (3-day period) will cause the 
hot Jupiter's orbit to precess completely in about 10 6 orbits 
= 3x 10 6 days, or about 8,000 years; a 3-Earth mass planet 
(m<i/M* ~ 10 -5 ) 1.5 times more distant than a hot Jupiter 
would cause precession in only about 10 5 orbits (800 years). 

The analysis above greatly underestimates the effect if 
resonant or near-resonant terms are important to the dy- 
namics. It is unclear whether resonant configurations are 
ubiquitous or happenstance; some extrasolar planetary sys- 
tems are already known to exhibit n ear-resonan t beha viour 
jRasio et al.lll992l; iFord et al.ll2005h . lAgol et al l |2005l ) and 
iHolman fc Murravl (|2005l )" include the effects of orbital reso- 
nances and the eccentricit y of the perturber's orbit; however, 
IHolman fc Murravl (2005) focus on the stochastic v ariation 
of the inter-transit interval and lAgol et al.l (|2005l ) do not 
consider secular terms in their calculations, which is the fo- 
cus of the analysis here. 

An explicit expression for the precession rate of the 
longitude of periastron, if the semimajor axis of the outer 
planet is much larger than that of the inner observed planet 
(a < a 2 ), is 



zu = 355 yr — — — — — 
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where P is the period of the inner planet. This asymptotic 
estimate of the rate of periastron advance (shown in Fig. [TJ 
significantly underestimates the rate for a > 1/4. In fact 
for a = 0.5 the estim ate falls short by a factor of two (c.f. 
lMiralda-Escudll2002l) . 




100 



Figure 1. The number of orbital periods of an inner planet re- 
quired for that planet's orbit to precess by 2n under the secular 
advance induced by an external perturber, times the mass ratio 
of the perturber to the central star (m,2/M t ). The a:-axis is the 
semimajor axis ratio of the outer to inner planet. For an unseen 
Jovian perturber the number of orbits for precession is the verti- 
cal axis time 10 3 and for a 3 Earth-mass perturber one multiplies 
by 10 5 . The dashed line is the approximation using the small-a 
limit of the Laplace coefficient. 



3 WHAT IS THE SENSITIVITY TO OTHER 
PLANETS IN THE SYSTEM? 

Fig. [2] shows the precession contribution from the three ef- 
fects on an observed transiting planet. The most interest- 
ing component is the contribution from other planets in the 
system, so we would like to have firm upper limits on the 
contributions from the other two effects in order to estimate 
the excess that might be due to an unseen planet. 

The relativistic contribution is the most certain. The 
mass of the star can generally be estimated to a few percent 
and the period of the orbit of the planet can be measured 
to a part in 10,000 or better; therefore, the value in Eq. [2] 
can be estimated to a high degree of precision. 

Since the GR-induced precession rate is simply propor- 
tional to the mass of the star, the fractional precision of this 
mass estimate will set the sensitivity limit of the periastron 
advance method. If we take 3% as a nominal M* precision 
then, if one observes a transiting hot- Jupiter with P—3 days, 
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1 10 100 

l/a = a s /a, 

Figure 2. The periastron precession rate of a hot Jupiter in a 
low-eccentricity, three-day orbit around a solar-mass star. The 
vertical dotted line shows the minimum separation of an exterior 
Jupiter- mass planet (see text). 



a precession rate different from the GR-rate by less than 3% 
will not provide a reliable detection of a planet or an oblate 
star; the 3-day period used for Fig. [2] gives by chance that 
the induced rate from the host star with J2 = 1CP 6 falls 
close to the detection limit. 

The value of the oblateness (J2) for the host stars of 
extrasolar planets is unknown - we don't even know it well 
for our Sun. However, we can estimate an upper limit to its 
value from the solar estimates and a scaling of the rotation 
rate of our sun to that of the host star if it is known. From 
observations of solar oscillations and theoretical considera- 
tions, the v alue of J2 for the sun is probably around a few 
times 10" 7 |Pireaux fc RozelotJl2003l ). IWinn et ail (|2005l ) re- 
cently determined the spin rate of the star in the transiting 
planetary system HD 209458. The value that they obtain 
tisini = (4.70 ±0.16) km/s is not much larger that the typi- 
cal values for the Sun of 1.4 — 2.0 km/s, so it is unlikely that 
the value of J2 for at least this system is much larger than 
that of the Sun. The value 10 -6 probably provides an up- 
per limit; this yields a oblateness contribution of about 3" 
per year, an order of magnitude less that the contribution 
from a nearby Earth-mass planet. It is unlikely that stars 
have values of J2 > 10~ 6 ; detected precession rates more 



than ~3" faster than the GR-induced rate would be strong 
evidence for the presence of an unseen planet. 

The detection of a periastron advance rate for an ex- 
trasolar planet with a three-day period that exceeded of the 
relativistic amount by 31" per year (Eq. (|3} would either 
indicate the presence of an unseen planet or a stellar J2 
of about 10 -5 , an order of magnitude larger than the up- 
per limits for the Sun. Either result is interesting. Assuming 
that the oblate star hypothesis is ruled out, such a detec- 
tion either implies (Fig. 01 a terrestrial-scale planet a few 
times further than the hot Jupiter, or a jovian-mass planet 
even more distant (for a given precession rate there is a 
one-parameter family of mass-distance for the external per- 
turber). 

The jovian-mass case would likely be limited by sev- 
eral constraints. Very large precession rates (many times the 
GR-induced precession) are unlikely to be foun d; if both 
planet s were jovian mass, then the analysis of iGladmanl 
1 19931 ) shows that an outer planet would not be stable if 
(<i2 — a\)/a\ < 0.24 (i.e., the outer planet's semimajor axis 
must be than >125% of the inner planet's or the system 
would be unstable; the vertical dotted line of Fig. [2}. Notice 
that a jovian perturber an order of magnitude further away 
than the hot Jupiter induces a precession rate comparable 
to the GR-induced rate. However, such a perturber should 
be trivially detectable already in the radial velocity data 
from the system; thus unless the host star is so noisy that 
radial-velocity techniques cannot be applied, the periastron 
technique is less sensitive to massive planets. 

However, while terrestrial-mass planets cannot be seen 
by the radial-velocity technique (the radial velocity signal 
is less than one meter per second for periods greater than 
three days), the periastron advance they would produce on 
a hot Jupiter is in principle detectable if they are up to ~ 5 
times more distant. 



4 HOW CAN ONE OBSERVE THE ADVANCE 
OF PERIASTRON? 

We can have the following information (loosely in order of 
difficulty) 

(i) Timing of the primary transit, 

(ii) Timing of the secondary transit and 

(iii) Timing of radial velocity data 



4.1 Primary transit 

As the orbit precesses the location of the transit relative to 
the orbit changes; specifically, the value of the true anomaly 
(v) at the center of the transit decreases at a rate of vj = 
8zu/P. To get the time of the transit, we have to relate the 
true anomaly v to the mean anomaly (A4 — the phase of the 
orbit increasing linearly in time from zero at one periastron 
to 27r at the next) through Kepler's equation, 

M=£-esin£, (13) 
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with 



tan- = 



1 — e v 

tan — 

1 + e 2 



(14) 



where £ is the eccentric anomaly. 

For simplicity we work in the frame rotating with the 
precessing orbit, so we take P to be the time between suc- 
cessive periastrons. The true anomaly that corresponds to 
the transit thus changes as the orbit precesses. The interval 
between two transits is simply the difference in the mean 
anomalies at each transit divided by the mean motion 2n/P. 
If the orbit is not precessing, this time is simply the period 
of the orbit. However, if the periastron advances, the timing 
between the transits is 



At = P 1 



where we have assumed that the periastron advance per or- 
bit is small. As the orbit precesses, the time between transits 
will change, according to 



dAt 
dt 
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At face value this gradual change of the timing of the tran- 
sits appears hopeless to detect because the change is propor- 
tional to the very small square of the advance of periastron 
per orbit. However, the difference in the timing accumulates 
from orbit to orbit as the orbit precesses, so in practice one 
would predict the timing of the transits from a few obser- 
vations and look for a difference from that prediction after 
many hundreds of orbits had passed. Essentially we are in- 
terested in the integral of At over many periods. 

If one observes several initial transits and determines a 
value of At, one can predict the timing of future transits. 
These predictions will be incorrect for a precessing orbit by 
the following amount 



tprcd - tactual = M(u) - M{vo) 



dM , A P 

AT ^-^) 2+ 
(18) 

where vo is the true anomaly during the transit at the initial 
epoch and v is the true anomaly during the transit at the 
later epoch. 

One can calculate 



d 2 M 
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or expanding for small eccentricities 
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From Eq. [17] one can see that for an error in the predictions 
to accumulate the second derivative the mean anomaly with 
respect to the true anomaly must not vanish. Thus, how 
quickly the error accumulates depends on the initial epoch 
of the observations. If the transit is initially occurring at 
periastron (y = 0) or apastron (y = 7r), the second derivative 
vanishes so it will take significantly longer for the time delay 
to become observable. 



To lowest order in the change in the periastron advance 
(or the true anomaly at transit), we have 
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so the delay accumulates quadratically in time. Using rea- 
sonable values for the various numbers we have 
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0.1 3 days \1 year 100" yr- 1 



(23) 



4-1.1 Error Analysis 



Eq. [23] makes it seem seem hopeless to detect the timing 
delay because one can determine the time of a particular 
transit to possibly ten seconds; therefore, naively one would 
expect to have to wait one hundred years before detecting an 
advance with t prc d — ^actual =10 seconds. Fortunately, one can 
detect the periastron advance in the series of transit times 
long before one could detect it in the timing of an individual 
transit. 

In practice one characterizes the timing of the transits 
with a formula of the following form 



t n = A+Bn+Cn 



P /c , 2 d 2 M 
t + At n+— (dzu) 

47r dv 2 



(24) 



where n is the number of the transit, to is the time of an ini- 
tial reference transit, Ato is the initial time between transits 
and the quadratic term contains the periastron advance. 
Using the standard results for \ 2 fitting, we obtain 
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To 
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(25) 
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where N is the number of the last transit sampled, ro is the 
fraction of transits with times and ao is the timing error on 
each transit. 

The upper limit obtained for the value of the advance 
per orbit Szu will be 



1/2 
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where we have ignored the fractional error in the values of 
P and d 2 M/du 2 . This yields 



1.6x10" 



N 
1000 



-5/4 



-1/4 



ao 3 day 0. 1 1 sin^ol 
Ws~ P e 



(Szu) GR — 5 x 10 



1 



1 



2/3 
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(28) 

as an upper limit on the advance per orbit. This is not much 
larger than the expected relativistic contribution of 



(29) 



so for TV greater than a few thousand the relativistic term 
will dominate over the statistical errors in the timing. These 
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estimates agree with the results of iMiralda-Escudi (|2002T l 
who considered the effects of periastron advance on the tim- 
ing of the primary transit and the duration of the primary 
transit. 

One thousand transits of a planet with a three-day orbit 
takes just a shade under eight and a quarter years. The up- 
per limit on dm decreases with time as t~ 5 ^ 4 until a reliable 
periastron advance is detected. After this time, the errors 
on this detection decrease as t~ s ^ 2 . The time to achieve the 
desired sensitivity scales as the P 3 / 5 so this technique is also 
applicable to planets with longer orbital periods. These es- 
timates assume that every transit is timed (ro = 1). The 
error analysis assumes that the observed transits are evenly 
spaced in time; it should be possible to devise an observing 
strategy that achieves errors similar to the ro = 1 case with 
many fewer observations — this is beyond the scope of this 
paper. 



where 

dM 
dv 



= nA 



: 1 — e cos £ 
1 + e cos v 



2e cos v + ^e 2 cos 2v 



1 — 2e cos v + ^e^ cos 2v — e A cos 3v + C(e 4 ) (32) 
If we take the difference between these two values we get 
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so the interval between two successive primary and two 
successive secondary transits differs by an amount propor- 
tional to the advance of periastron per orbit. This should be 
compared with observations of the primary transit alone in 
which the advance of periastron only enters at second order. 

Furthermore, the difference in the time between the pri- 
mary and secondary transits also contains some valuable in- 
formation. We have 



4-1-2 Why don't observations of the primary transit tell 



tjT 



us morei 



Each time the planet orbits the star is takes a bit less than 
an orbital period for it to reach the point of primary transit, 
because the orbit is shifting a bit. However, since we don't 
know the radial orbital period itself, this time is essentially 
unobserved. The time for the planet to cover the missing an- 
gular distance is related to the distance from the star to the 
planet at transit. As the orbit precesses, this distance will 
change which in turn will change the time between transits. 
It is this change in the time between transits that we try 
to observe. The correction in the time between transits is 
proportional to the periastron advance. The change in the 
distance between the star and the planet from orbit to orbit 
is also proportional to the periastron advance. Combining 
these facts indicates that the change in the time between 
transits is second-order in the small periastron advance; con- 
sequently, it takes a relatively many orbits to detect the pe- 
riastron shift, if one times only one type of transit. 



4.2 Secondary transit 

Looking at the secondary transit (when the planet goes be- 
hind the star) does not just provide a new set of times to fit 
but also provides new information and possibly a faster way 
of detecting unseen planets in the system. The secondary 
transit occurs when the true anomaly is 180 degrees away 
from where the primary transit occurs; therefore, we will de- 
note quantities that describe the secondary transit with the 
subscript w . The primary transit is given by the subscript o- 
Let us examine at the time between two successive primary 
and two successive secondary transits. Using the earlier for- 
mulae we have to lowest order in the advance of periastron 
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If the orbit is eccentric this will differ from half of the orbital 
period. Expanding in the eccentricity we have 



in — to 
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The first term in the series is twice the value of the first 
term in the series for d 2 M/dv 2 , so the interval between the 
primary and secondary transits helps to calculate the pe- 
riastron advance when one uses the timing of the primary 
transits. 



4-2.1 Error Analysis 

How well can we determine the time of the transits and the 
time between successive transits? Fitting the transit times 
to a timing model 



and 



to,n = to + Aton 



£"7T,n — t-n ~\~ ^i-ftTh- 



(36) 



(37) 



Because the periastron advance now enters in the difference 
between the interval between the successive transits, we only 
need to fit the times to first order in the number of the 
transit "n". From the x 2 — analysis we obtain the following 
error estimates 



and 



2j2(ryJV) + 1] 
(r*N)[(r*N) - 1] 



12< 



1/2 



1/2 



(38) 



(39) 



where rv and ov are the fraction of secondary transits with 
times and the error in the timing of the secondary transit. 
The error in the time of the initial transit scales as iV~ 1//2 
where TV is the number of orbits that have elapsed between 
the first and last one observed. 
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We are interested in the differences 

p 

Ato-AU = (-4ecos^ - 2e 3 cos3j/ + 0(e 5 )) 8vj— (40) 

Z7T 



and 



<XAt -At„ 



12 



1/2 



(41) 



for large N, and especially the error in 



17 6 z 



7 x 10" 



0.1 / TV 



(10 



e| cos i/o | V 1000 
a 2 . 

+ 



-3/2 



3 day 



(10 s) 



1/2 



(42) 



Combining results for the secondary transit with those from 
the primary transit yields an increase in sensitivity of a fac- 
tor of ~ 200. If we could time the secondary transit to the 
same precision of ten seconds it would take only N ~ 100 
to detect an Earth-like planet within twice the semimajor 
axis of the observed planet. With current instruments and 
the brightest targets, the secondary transits can be timed 
to a precision of about 100 s (J. Matthews, priv. comm.), 
yielding an estimate of about N ~ 400 orbits for a sim- 
ilar detection. The time to achieve the desired sensitivity 
scales as P 1/3 so this technique is also applicable to planets 
with longer orbital periods; furthermore, one is sensitive to 
smaller planets in systems with larger eccentricities. 



4-2.2 Why does the secondary transit help so much? 

After analyzing the primary transit, one saw how difficult 
it was to disentangle the change in the angle of periastron 
from the observations of the orbital period of the system. 
The timing of the secondary transit breaks this degeneracy, 
and it is straightforward to understand why. Unless the or- 
bit is perfectly circular (or if we are so unlucky as to have 
vq w 7t/2 or 37r/2 in the epoch of observations), the planet 
is at different distance from the star at the primary and sec- 
ondary transit, so according to Kepler's Second Law (con- 
servation of angular momentum) its angular velocity along 
the orbit is different at these two times. If the periastron 
shifts as the planet orbits, it takes a different amount of 
time to cover the missing angle at the primary than at the 
secondary transit; consequently, the time between secondary 
transits differs from that between the primary transits — if 
one can detect this time difference one can detect the ad- 
vance of periastron. 



4.3 Radial Velocity Information 

The radial velocity information is arguably the most difficult 
to obtain. It turns out that it is essentially the least useful 
(at least in quantity) for the purposes of characterizing the 
periastron shift. It is difficult to imagine obtaining timing of 
the radial velocity data with a precision of tens of seconds, so 
it is not directly useful in getting additional timing points, as 
we did with the secondary transit. In principle, one would 
find that the time interval between when the star passed 



through a particular radial velocity and when it repeated 
itself would depend on the radial velocity in question. 

However, the period found by fitting the radial velocity 
curve is typically precise to about one second . This time 
interval would only differ from the interval between tran- 
sits by a tiny amount, on the order of the periastron ad- 
vance. Determining accurately the relationship between the 
time between periastrons (what we have called the period) 
and the period found by fitting the radial velocities requires 
Monte Carlo simulations of the observed data. One can also 
gain some insight into what time interval emerges from fit- 
ting radial velocities by considering orbits that are nearly 
circular. 

When one fits the radial velocity measurements one is 
most sensitive to parts of the orbit with large accelerations 
to or from the observer. The acceleration along the line of 
sight reaches an extreme when the jerk vanishes. To first 
order in the orbital eccentricity 



d?V\ os 

dt 2 



sin (y — !/()) + 



2e (11 cos v sin (y — vq) + 5 sin z/q) 



so the jerk vanishes where 



sin (y — vo) = — lOe sin vq + 0(e 



(43) 



(44) 



To lowest order in the eccentricity, the radial velocity mea- 
surements are equally sensitive to the timing at the primary 
(y — Uq = 0) and secondary transits (y — vo = 7r), so we as- 
sume that the time interval determined by fitting the radial 
velocity is given by the average of the two intervals discussed 
earlier 



At RV = i (At + At„ 



+ 0(e ) =P 1 



(45) 

The "period" obtained by fitting the radial velocity data dif- 
fers slightly by the period between periastrons or the period 
between primary transits. We confirmed this by generating 
radial velocity data with an advancing periastron and fit- 
ting these data with purely Keplerian radial velocity curves. 
These simulations gave Eq. US] for small eccentricities. 

If we take the difference between the two observable 
quantities we get 



Au 



Ato = 2e cos ;/o Svj — . 

2n 



The error in this quantity is given by 



O "a r 



°"At RV +°0 



12 



1/2 



2tt 



roiV 3 / P 2e|cosi/o| 



10" 



_2 

(100 ms) 2 

3 day 0.1 
P e\cosv \ 



+ 0.12r " 1 



(10 s) 



(— 

Vioo 



(46) 



(47) 



1/2 



(48) 



We see that the timing errors in the radial velocity mea- 
surements dominate over the transit timing for a quoted 
precision in the radial-velocity period of 100ms. Currently 
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the best period estimates for hot Jupiters without t ransit 
information are precise to 800ms (|Butler et all 12006), but 
analysis of radial velocity measurements over longer base- 
lines would provide a more precise estimate of this period. 

Even without a detailed understanding of the relation- 
ship between the velocity and transit timing, the radial ve- 
locity data is crucial to convert a observed timing solution 
into a periastron shift by determining the values of the ec- 
centricity and the true anomaly at transit, and in combi- 
nation with the timing of either the primary or secondary 
transit could yield hints of the periastron advance due to 
other planets in the system. 

As the orbit precesses, the radial velocities observed 
over a orbit will also shift. However, over the time required 
to detect a Earth-like planet the orbit will only precess about 
an arcminute. It is difficult to imagine that radial velocity 
measurements will become so sensitive as to characterize 
an orbit to the required precision of 10 -4 . If they did, one 
could probably detect the planet causing the precession in 
the radial velocity data already. 



4.4 Rapid Precession 

If the transiting planet's orbital eccentricity goes to zero, 
the precession rate of its pericenter longitude will formally 
go to infinity, and Eq. [23] misleadingly indicates that the 
pericenter rate will be trivial to detect. However, Eg. 1151 is 
not correct in the case of rapid precession. 

If the rate of the precession is a constant (Suj) per radial 
orbit we have following equation for the true anomalies of 
the transits 



limited by e 2 P/(47r). We have 



v mod 2-7T 



v 



mod 2-7T. 



(49) 



If we assume that the orbit is a precessing ellipse and that 
the angular momentum of the observed planet is conserved, 
we can use Kepler's equations to determine the time cor- 
responding to each true anomaly. To lowest order in the 
eccentricity we have 



,v_ _ P_ £ 

2^ ~~ 2^T 



sin 2v 



(50) 



If we first ignore the eccentricity, we find that the time be- 
tween two successive transits is given by the angular period 



At m = P ( 1 



and 



Szu 
~2t7 



Av = v 2 — V\ = —duj [ 1 + — 



(51) 



(52) 



Looking at Kepler's equation we find that the correction to 
this quantity introduced by the eccentricity of the orbit is 



At<°> 



2tt 4 



(sin 2i/2 — sin 2vi) 



(53) 



= At 



1 + — ^ — [sin Av cos (2vi + Av)] 1 (54) 

2-7T J 4-7T J 

For a perturber in an elliptical orbit Szu is inversely 
proportional to the ecce ntricity of the observed planet 
|Murrav fc Dermottl l200(il ) ; thus Eq. [54] indicates that the 
observable correction to the transit time is proportional to 
the eccentricity, proving that the timing error becomes un- 
observable as e tends to zero rather than diverging. 



5 SPECIAL RELATIVISTIC CORRECTIONS 

The foregoing analysis focused on the angles necessary for a 
transit to occur. It was essentially geometry with any kine- 
matics. Specifically it neglected the time for light to travel 
across the system. The variation in the distance of the planet 
and the star from transit to transit as the orbit precesses 
would affect the times that we observe the transits to occur. 

The light travel time will affect the observed time differ- 
ence between the primary tr ansit and t he secondary transit 
that immediately follows it l|Loebll2005l ). 



[*» - ioloba = [M{y + tt) -M(v )] — + 



2?r 



where 



- [r(yo) + r(v + tt)] 
c 



a(l 



1 + e cos v 



(55) 



(56) 



The light-travel time will cause the primary transit to ap- 

earlier. The 



pear to occur about 20 s ^; 



2/3 



\ M eJ 



1/3 



3 days 

appearance of the secondary transit will be delayed by simi- 
lar interval. This timing signature may be used to constrain 
the physical size of the orbit; however, the difference in the 
inter-transit interval due to geometry is on order of the pe- 
riod of the orbit, a factor of 10 4 larger, so the eccentricity 
of the orbit must be known accurately for the relativistic 
corrections to be useful. 

Other quantities that we have examined are the the in- 
terval from primary to primary transit and from secondary 
to secondary transit, and the time derivative of these quan- 
tities. If the orbit did not precess the light travel time would 
not affect either of these intervals, so we know that the rel- 
ativistic correction to these intervals will be proportional to 
the change in a angle of periastron during an orbit. 

We have 



and 



[A*o] ob 



[AM ob 



1 



P_ dM 

2tt dv 

P_ dM 

2n dv 



1 dr 

c dv 



+ 



1 dr 
dv 



(57) 



(58) 
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The ratio of the two corrections to the interval between 
primary transits is 



tan vq a 



At: 



0,ang 



8 x 10 



cP 

_ 5 tan vq ( M» 3 days 

TV 1 M© P 



1/3 



(59) 



(60) 



The rati o is similar f or the secondary transits; therefore, as 
found bv lAgol et ail (|2005l ) it is safe to ignore the relativistic 
corrections to the intervals between two similar transits. 



6 OUTLOOK 

6.1 Systems with a transiting planet 

Nearly two hundre d extrasolar planets have been detected as 
of the beginning of l2006l l|Butler et al . 20060 The eccentric- 
ity and the longitude of periastron have been measured for 
most of these planets; unfortunately, for the few transiting 
planets, this information is lacking, so we will use the orbital 
elements for all the nearby exoplanets to calculate the sen- 
sitivity of timing measurements to periastron advance and 
more importantly how long of a timing series would be re- 
quired to achieve a precision of 10 -7 in the measurement of 
S-oo; this is sufficient to detect an Earth-mass planet orbiting 
within a factor of four of the semimajor axis of the transiting 
planet for a one-solar- mass host star. Fig. [3] shows that for 
about ten systems this level of sensitivity could be achieved 
within a decade (these results have assumed that the sam- 
pling rate r is unity; determining the optimal sampling rate 
and schedule considering observational constraints is beyond 
the scope of this paper). 

To be certain of the presence of an unseen planet one 
has to be sure that the periastron precession exceeds that 
expected from the star. For this discussion we shall assume 
that the contribution due to stellar oblateness is small. Ac- 
cording to Fig. [2] even for a three-day period, the expected 
contribution from oblateness is a factor of 30—300 below 
the relativistic value. This ratio increases as P 2 / 3 so oblate- 
ness does not make the dominant contribution for any of 
the systems in Fig. [3] On the other hand, Fig. [4] shows the 
expected relativistic precession is greater that 10 -7 for most 
of the systems in the catalogue; consequently, the mass of 
most of the host stars for these systems must be determined 
to better than about a tenth of a solar mass. Otherwise, the 
error in the determination of the mass of the host star will 
dominate the statistical error in the timing. 

Fig. [5]shows the sensitivity to finding additional planets 
in the transiting systems with a ten-year time series assum- 
ing that the mass of the host star is known with a precision 
of 0.03Mq. The properties of the t ransit itself , especially 
the shap e of the ingress a nd egress l|Loebll2005h and dura- 
tion fe.g lWinn et al.ll2006l ). may constrain the mass of the 
star further. Four objects stand out in Figs. [3HS1 GJ 436 b, 
HD 118203 b, HD 33283 b and HD 74156 b require the 



1 We used the updated catalogue at http://exoplanets.org 
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Figure 3. The duration of the time series of primary and sec- 
ondary transits (or primary transits and radial velocity timing) 
to yield a detectable 5m of 10 — 7 for the planets in the catalogue 
of lButler etHI {2006), assuming a ten-second error in the timing. 
This is the sensitivity required to detect an Earth-mass planet at 
the three-sigma level within a factor of four in semimajor axis 
of the transiting planet, orbiting a solar mass star. The required 



length of the time series increases as a 



-2/3 



shortest time series to achieve a sensitivity of 10 -7 in 5w. 
These systems have moderate eccentricities between 0.2 and 
0.6. Because of its short 2.6-day orbit, the orbital preces- 
sion of GJ 436 b would be dominated by relativistic effects 
with (<5-to )gr ~ 2.3 x 10~ 6 assuming a mass of O.41M0 for 
GJ 436 (|Butler et al.ll2006l ). There are several planets in 
the catalogue whose orbits will precess at a rate similar to 
GJ 436b due to GR, but the geometry and eccentricity of 
the orbits are not as favourable for detecting the preces- 
sion by transit timing. On the other hand, the maximum 
sensitivity to unseen planets are found in the HD 74156 
and HD 168443 systems because their relatively long orbits 
(52 days and 58 days) reduces the expected contribution 
due to GR and stellar oblateness. The minimum detectable 
value of (<5ti7)pi anc t is about 5 x 10~ 8 , sufficient to detect an 
Earth-mass planet with a period less than about 460 days 
after a decade of observation. 

Unfortunately, those planets that require the least time 
to constrain additional bodies in the system have relatively 
long periods, on the order of 50—100 days. It is a factor of 
seven to ten times less likely that planets with such long pe- 
riods transit their star than the hot Jupiters with three-day 
periods. However, as more transiting planets are discovered 
one would expect to find both planets with longer periods 
and higher eccentri cities than th e curre nt cohort of transit- 
ing planets. Indeed [Butler et alj |2006) list several planets 
that orbit their stars with periods of a few days in eccen- 
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Figure 4. The advance of th e long itude of periastron per orbit 
for the planets of lButler et al,l d2006l ) from general relativity. The 
four systems for which the relativistic advance is most rapid are 
HD 86081, HD 149026, tau Boo, HD 73256. 



trie orbits (e > 0.1). None of these systems have exhibited 
transits, but it is only a matter of time before a transiting 
planet with an eccentric orbit is found and these techniques 
may be brought to bear. 

6.2 Systems with more than one observed planet 

A system that already has more than one planet detected 
provides a chance to characterize an additional planet and 
verify that the stellar oblateness does not dominate over 
the planetary signal. Of course, each observed planet will 
induce a periastron precession in the other, thus giving the 
planetary masses directly. This must be subtracted from the 
observed shifts along with the relativistic shifts. Because 
both planets orbit the same star the contribution due to 
the stellar oblateness is proportional P~ 7 ^ 3 ; consequently, 
unless the residual periastron precession rate in each planet 
that remains is proportional to P _7//3 , there must be other 
planets in the system. 

If one can argue from other data that the value of J2 
is small, one can use the residual precession rate to find the 
mass of the unseen planet and its semimajor axis. With a 
single observed planet one can only constrain the combina- 
tion the determines Suj. With three observed planets, one 
could unravel J2 and the properties of an unseen planet or 
alternatively the properties of an unseen planet and the pos- 
sibility of further planets! 




1 10 100 1000 io- 1 10 5 

Period [days] 



Figure 5. The sensitivity to the advance of the longitude of peri- 
astron caused by ad ditional plane t s (or stellar oblateness) in the 
planetary systems of lButler et al . (2006) assuming that the mass 
of the host star is known to O.O3M0 with a ten-year baseline of 
transit timing. To detect an Earth-mass planet with an orbit up 
to a factor of four larger than the observed planet the sensitivity 
to the planetary contribution to the precession must be less than 
IO" 7 . 



6.3 Systems with fewer than one observed planet 

It may seem odd to suggest measuring the periastron ad- 
vance in systems without any planets yet discovered. How- 
ever, the techniques outlined here could be used in eclipsing 
binary systems. In single- lined systems it could be used to 
constrain the total mass of the system through the rela- 
tivistic term. In double-lined systems or if the total mass 
of the system is known, measurements of the periastron ad- 
vance through timing of the radial velocity data and one or 
both transits could be used to constrain the values of J2 for 
the stars or to look for planets in the systems. One could 
search for planets with the long-term timing measurements 
of eclipsing binaries that have already been taken. 



7 CONCLUSIONS 

Accurate timing data (of pulsars) allow ed the discovery of 
the f irst extrasolar Earth-mass planets (|Wolszczan fc Frail 
1992). Careful and accurate timing of planet transits and 
radial velocity data is sensitive to additional planets down 
to the mass of Earth and below. The combination of two 
sets of timing data provides much stronger constraints on 
the presence of additional bodies in the system than look- 
ing at the primary transits alone (c.f. lMiralda-Escude1l2002l ; 
ISchneiderll2uO^ ; lAgol et al.|[2005l ; iHolman fc Murray! 120051 ) . 
The timing signature of an Earth-mass planet is an induced 
shift in the periastron of the orbits of the known planets. 
There are generally two dominant contributions to this shift: 
general relativity and other planets. The stellar oblateness 
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can also contribute but only competes with the other effects 
if the oblateness is nearly two orders of magnitude larger 
than that of the Sun. Consequently, if the observed peri- 
astron shift exceeds the relativistic expectation, either the 
system has additional planets or the parent star has an un- 
usually large oblateness. Even the less likely possibility of a 
large oblateness would give tantalizing hints to the origins 
of these close-in extrasolar planets. More likely would be the 
presence of a Earth-mass or even a Mars-mass planet in a 
nearby orbit to the known planet. 
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